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Abstract.  We  propose  a  rate  process  model  for  describing 
photochemical  damage  to  retinal  cells  by  short  wavelength 
laser  exposures.  The  rate  equation  for  photochemical  dam¬ 
age  contains  a  positive  rate  that  is  temperature  indepen¬ 
dent,  and  a  negative  (quenching)  rate  that  is  temperature 
dependent.  Using  the  traditional  Arrhenius  integral  to  de¬ 
scribe  thermal  damage,  we  derive  damage  threshold  doses 
for  both  thermal  and  photochemical  mechanisms,  and  show 
that  the  model  accounts  for  the  sharp  transition  from  ther¬ 
mal  to  photochemical  damage  thresholds  that  have  recently 
been  observed  in  an  in-vitro  retinal  model.  © 2011  Society  of 
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1  Introduction 

For  long  pulse  durations  (>1  s),  laser  induced  tissue  damage 
can  be  divided  into  two  broad  categories:  thermal  and  nonther- 
mal  (photochemical).  Thermal  damage  has  been  modeled  by 
the  Arrhenius  damage  integral,  first  proposed  by  Henriques1  in 
1947,  with  great  success.2-5  Other  models  have  been  developed 
to  address  the  shortcomings  of  the  Arrhenius  integral,6,7  but  it 
continues  to  be  a  popular  method  for  predicting  thermal  damage 
thresholds  due  to  its  relative  simplicity. 

No  such  damage  model  exists  for  predicting  photochemi¬ 
cal  damage  thresholds.  Perhaps  the  demand  for  such  a  model 
has  been  low  because  photochemical  damage  is  defined  by  the 
principle  of  reciprocity,  where  the  total  radiant  exposure  re¬ 
mains  constant  with  exposure  duration,  allowing  damage  to  be 
predicted  by  simply  considering  the  total  energy  delivered. s 
However,  in  his  in-vitro  retinal  model,  Denton  et  al.  showed 
that  at  413  nm,  damage  thresholds  follow  a  thermal  trend  line 
for  exposure  durations  up  to  60  s  and  can  be  a  factor  of  2 
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higher  than  the  reciprocity  radiant  exposure.4  From  this  work, 
the  damage  threshold  mechanism  seems  to  transition  sharply 
from  thermal  to  photochemical  somewhere  between  60  and 
100  s,  and  a  similar  trend  has  been  described  in  an  animal  model 
as  well.10 

In  this  paper,  we  propose  a  rate  equation  to  describe  pho¬ 
tochemical  damage  that  includes  a  temperature  dependence. 
The  model  assumes  the  creation  of  some  activated  molecule 
(B*)  through  a  photon  absorption  that  is  independent  of  tem¬ 
perature  and  a  temperature  dependent  quenching  rate  that  re¬ 
moves  B*  (Fig.  1).  We  treat  thermal  and  photochemical  damage 
mechanisms  as  independent  processes,  the  damage  threshold 
is  defined  as  the  lowest  threshold  of  the  two.  The  temperature 
independence  of  B*  creation  and  the  temperature  dependence 
of  quenching  leads  to  a  sharp  cutoff  for  photochemical  damage. 
Above  some  critical  temperature,  the  photochemical  damage 
mechanism  is  completely  negated,  allowing  the  damage  thresh¬ 
old  to  be  determined  by  thermal  damage  at  energies  much  higher 
than  the  reciprocity  energy. 


2  Damage  Model 

2.1  Thermal  Damage  Rate 

We  assume  thermal  damage  to  be  described  by  the  usual  Arrhe¬ 
nius  rate  equation,  where  some  molecule  A  (typically  assumed 
to  be  a  protein)  is  converted  to  a  denatured  state  A* .  The  incident 
light  provides  a  heat  source  which  raises  the  temperature  of  the 
cells,  increasing  the  rate  of  this  reaction.  The  damage  rate  is 
given  as 


dt 


(i) 


Here  the  brackets,  [•],  denote  molar  concentrations,  and  R  and 
T  (t)  are  the  gas  constant  and  system  temperature  (in  Kelvin), 
respectively.  The  two  constants,  Eai  and  Z\,  are  termed  acti¬ 
vation  energy  and  collision  frequency.  These  rate  coefficients 
must  be  experimentally  determined,  and  are  strongly  dependent 
on  the  system  and  damage  endpoint  assessed.  Several  sets  of 
coefficients  exist  in  the  literature.11 


2.2  Photochemical  Damage  Rate 

We  propose  a  two-process  rate  equation  to  describe  photochem¬ 
ical  damage.  Some  molecule  B  absorbs  a  single  photon  to  create 
an  activated  molecule  B*.  This  activated  molecule  may  then  be 
deactivated  through  an  interaction  with  some  other  molecule, 
§,  provided  by  the  cell.  The  activated  molecule,  B*,  may  also 
initiate  a  chain  of  events  that  eventually  leads  to  cell  death 
(Fig.  1). 

The  total  rate  of  creation  for  B*  can  then  be  written  as  a 
positive  rate  dependent  only  on  the  incident  photon  flux  plus 
a  negative  rate  that  accounts  for  the  chemical  reaction-based 
removal  of  B*, 

=  e<h[S]  -  Z2e-,^/RT(l,\^\\B*\.  (2) 

dt 

Here,  4>  is  the  incident  photon  flux  density  (number  of  photons 
per  area  per  time),  e  denotes  an  efficiency  factor,  which  is  the 
probability  that  an  incident  photon  is  absorbed  and  creates  the 
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Fig.  1  Photochemical  damage  mechanism.  A  single  photon  (y)  is  ab¬ 
sorbed  by  B  to  produce  B*  and  some  other  molecule  £  reacts  with  B* 
to  regenerate  B,  or  a  non-reactive  product  of  B.  Any  accumulation  of 
B*  can  be  available  to  trigger  a  chain  reaction  resulting  in  cell  death 
(dashed  line) 


activated  molecule,  B* ,  and  Ea 2  and  Z2  are  activation  energy 
and  collision  frequency  associated  with  £  and  B*  reacting  to 
give  B.  We  will  assume  that  the  molecules  B  and  £  are  abundant, 
so  their  concentrations  can  be  considered  constant.  We  then  use 
the  simplified  rate  equation  for  B* 

1  =  e<J>  -  Z2e~Ea2,RT<,)[B*].  (3) 

dt 

We  assume  photochemical  damage  to  be  caused  by  a  buildup  of 
B*,  and  so  the  damage  threshold  will  correlate  to  a  concentration 
of  B*.  The  relationships  between  positive  and  negative  reaction 
rates,  the  temperature  independence  of  damage  creation,  and  the 
temperature  dependence  of  damage  quenching,  are  unique  to  the 
photochemical  mechanism. 


3  Damage  Threshold  Predictions 

We  seek  to  find  an  expression  for  the  threshold  radiant  exposure 
for  causing  damage  to  cells.  In  our  analysis  here,  the  thermal  and 
photochemical  damage  mechanisms  compete.  Thermal  dam¬ 
age  is  caused  when  A*  reaches  some  threshold  concentration, 
[A*](h,  and  photochemical  damage  is  caused  when  B*  reaches  a 
threshold  concentration,  [B*]th.  Both  thermal  and  photochemi¬ 
cal  mechanisms  will  have  a  threshold  radiant  exposure  to  cause 
damage.  The  two  rate  equations,  (1)  and  (3),  are  integrated  to 
give  [A*]  and  [B*]  as  a  function  of  time  and  then  solved  for  the 
radiant  exposure  (H)  that  gives  [A*]th  or  [B*]th  for  the  expo¬ 
sure  duration  r.  The  lower  of  these  two  thresholds  will  be  the 
observed  damage  threshold  for  the  exposure. 

In  the  experiments  conducted  by  Denton  et  al., 9  retinal  pig¬ 
ment  epithelial  (RPE)  cells  were  exposed  to  413  nm  laser  radi¬ 
ation  for  durations  of  0.1,  1.0,  20,  40,  60,  100,  and  200  s.  For 
exposures  greater  than  about  20  s,  the  system  reached  steady- 
state  early  in  the  exposure,  and  the  temperature  can  be  assumed 
constant  throughout.  The  transition  from  thermal  to  photochem¬ 
ical  damage  occurred  somewhere  between  60  and  100  s,  which 
we  seek  to  describe,  so  we  can  assume  the  exposure  temperature 
to  be  constant. 

Due  to  the  linearity  of  the  heat  equation,  the  peak  steady- 
state  temperature  can  be  written  as  a  linear  function  of  the  peak 
incident  power,  E0, 


T  =  mE0  +  T0, 


(4) 


3.1  Thermal  Damage  Threshold 

Integrating  Eq.  (1)  with  the  temperature  held  constant  gives 

[A*]  =  Zie-E°'/RT  r,  (5) 


where  r  is  the  exposure  duration  and  T  is  the  exposure  tem¬ 
perature.  The  threshold  temperature,  7*,  is  the  temperature  that 
leads  to  the  threshold  damage  accumulation,  [A*]^.  However, 
the  exposure  temperature  is  related  to  the  exposure  irradiance 
through  Eq.  (4).  Inserting  T  =  mEt j,  +  Tq  into  Eq.  (5)  and  solv¬ 
ing  for  E,i,  gives  the  damage  threshold  irradiance.  Most  rate 
coefficients  (Zj  and  Ea\  )  for  thermal  damage  are  measured 
assuming  that  [A*]th  =  1,  so  the  damage  threshold  irradiance  is 


f  Ah  —  — 

m 


Ea\ 

Bln(Z,r) 


(6) 


3.2  Photochemical  Damage  Threshold 

To  integrate  Eq.  (3),  we  define  K2  (T)  =  Z2e~Eal^RT  and  mul¬ 
tiply  by  the  integration  factor  eKlt  which  gives, 

[B*]=^[  (7) 

We  must  define  [B*]*  to  obtain  a  damage  threshold  exposure. 
For  this,  we  turn  to  the  principle  of  reciprocity.  Here,  the  radi¬ 
ant  exposure  remains  constant  as  the  exposure  duration  varies. 
It  follows  that  as  the  exposure  duration  increases  the  incident 
irradiance  is  reduced,  which  in  turn  implies  that  the  exposure 
temperature  is  lowered.  Reciprocity  suggests  that  the  quenching 
rate,  K2  (T ),  is  small  for  low  temperature  rises. 

We  define  the  reciprocity  radiant  exposure  as  H, .  The  photon 
flux  density  can  be  written  in  terms  of  the  radiant  exposure  as 
<t>  =  H/rhv,  where  hv  is  the  photon  energy  and  h  is  Planck’s 
constant.  At  reciprocity  then,  Eq.  (7)  can  be  written  as. 


[B*k 


6  Hr 
K2  ( T )  r Irv 


(8) 


Since  K2  (T)  is  small  here,  we  may  expand  the  exponential  in 
a  power  series  and  ignore  all  second  order  and  higher  terms 
{(K2t)2  1)  and  Eq.  (8)  simplifies  to 

*  cHr 

[B  ].h  =  (9) 

hv 


Using  this  definition,  we  can  solve  Eq.  (7)  for  the  threshold 
radiant  exposure, 


Hth  =  HrK2(T)r[l-e-K^T)t]  '  .  (10) 


This  gives  the  damage  threshold  radiant  exposure  for  some  fixed 
temperature.  However,  the  exposure  temperature  is  also  depen¬ 
dent  on  H  through  Eq.  (4)  so  Eq.  (10)  becomes 

Hth  =  H,  K2  (m^  +  To)  t  [1  -  e-*4m/Wr+r„]r]-' 

(11) 

We  use  a  numerical  search  algorithm  to  determine  the  radiant 
exposure,  //,h.  that  satisfies  this  equation. 


where  7o  is  the  initial  (ambient  or  body)  temperature.  The  slope 
of  this  line,  m,  will  depend  on  the  specific  laser  parameters  and 
system  configuration. 


4  Results 

Measurements  for  Z\  and  Ea\  exist  in  the  literature,  but  no  such 
measurements  exist  for  Z2  and  Ea2.  However,  we  can  use  the 
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Fig.  2  Predicted  temporal  action  profile  (line)  compared  to  the  mea¬ 
sured  damage  thresholds  by  Denton  et  al.'J  (points).  The  inset  shows 
photochemical  damage  predictions  for  three  different  collision  fre¬ 
quencies:  1080  (top),  Id100  (middle),  and  1 01 20  (bottom). 


deliver  the  same  number  of  photons.  But  the  increase  in  irradi- 
ance  leads  to  a  larger  temperature  rise,  which  eventually  disables 
the  photochemical  damage.  We  have  not  attempted  to  identify 
the  underlying  pathways  responsible  for  causing  photochemical 
damage;  the  model  only  assumes  that  the  damage  rate  can  be 
described  by  a  rate  equation.  A  discussion  of  possible  photo¬ 
chemical  mechanisms  at  the  molecular  level  can  be  found  in 
Ref.  13. 

The  damage  threshold  predictions  reproduce  the  sharp  tran¬ 
sition  between  thermal  and  photochemical  damage  observed  in 
the  in-vitro  data  reported  by  Denton  et  al.9  The  two  rate  coef¬ 
ficients  for  photochemical  damage  quenching  have  been  esti¬ 
mated,  but  these  would  need  to  be  determined  experimentally 
to  make  quantitative  predictions.  In  principle,  these  coefficients 
would  need  to  be  determined  for  each  tissue  and  endpoint  to 
be  described,  but  could  then  be  used  to  predict  whether  or  not 
photochemical  damage  is  caused  under  various  exposure  condi¬ 
tions,  similar  to  how  the  thermal  damage  Arrhenius  coefficients 
are  used  today. 


critical  temperature  (the  temperature  for  which  the  Arrhenius 
rate  is  equal  to  one11),  Tc  =  Ea2/ R  In  Zi,  to  define  the  relation¬ 
ship  between  Z2  and  Ea 2.  We  previously  showed11  that  the  ex¬ 
pected  temperature  for  the  100  s  threshold  exposure  was  around 
38°C,  while  the  expected  temperature  for  the  60  s  threshold  ex¬ 
posure  was  about  48°C.  This  is  in-line  with  the  rule  of  thumb 
that  thermal  damage  requires  about  a  10°C  temperature  rise. 
Here,  we  assume  that  the  critical  temperature  associated  with 
the  photochemical  rate  process  is  47°C  (320.15  K). 

We  use  the  thermal  damage  rate  coefficients  reported  by 
Welch  and  Polhamus  for  thermal  damage:12  Ea\  =  6.28 x  105J 
mol-1,  Z\  =  3.1xl0"s-1.  The  expected  temperature  rises 
reported  by  Denton  et  al.9  give  m  =  0.26,  and  the  reciprocity  en¬ 
ergy  was  940  J  cm-2.  We  have  computed  damage  threshold  pre¬ 
dictions  for  three  different  values  of  Z2:  1080,  1010°,  and  10l2°. 
With  Ea 2  =  320. 15R  In  Z2,  this  gives  4.9x  105,  6.13x  105,  and 
7.35  x10s  for  Ea2-  The  model  predicts  the  same  sharp  transi¬ 
tion  from  thermal  to  photochemical  damage  as  was  observed 
(Fig.  2),  but  the  transition  region  depends  on  the  value  of  Zi. 
The  smallest  collision  frequency  (Z2  =  1080)  gives  photochem¬ 
ical  damage  at  exposure  duration  less  than  60  s,  while  the  largest 
(Z2  =  10120)  leads  to  no  photochemical  damage  at  100  s.  The 
best  agreement  was  with  Z2  =  101(,(). 


5  Conclusions 

We  have  presented  a  theoretical  model  that  assumes  a  temper¬ 
ature  dependent  photochemical  damage  rate  and  provides  an 
explanation  for  the  sharp  transition  between  thermal  and  pho¬ 
tochemical  damage  recently  observed  in  in-vitro  and  in-vivo 
data  for  short  wavelength  damage  thresholds.  As  the  expo¬ 
sure  duration  is  decreased,  the  irradiance  must  be  increased  to 
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